Ising model in scale-free networks: A Monte Carlo simulation 



O 

o 

(N 
o 

Q 



Carlos P. Herrero 

Instituto de Ciencia de Materiales, Consejo Superior de InvesUgaciones Cienttficas (CSIC), Campus de Cantohlanco, 

Madrid, Spain 
(February 2, 2008) 

The Ising model in uncorrelated scale-free networks has been studied by means of Monte Carlo 
simulations. These networks are characterized by a degree (or connectivity) distribution P{k) ~ k 
The ferromagnetic-paramagnetic transition temperature has been studied as a function of the pa- 
rameter 7. For 7 > 3 our results agree with earlier analytical calculations, which found a phase tran- 
sition at a temperature T'c(7) in the thermodynamic limit. For 7 < 3, a ferromagnetic-paramagnetic 
crossover occurs at a size-dependent temperature Tco, and the system is in the ordered ferromagnetic 
state at any temperature for a system size — > 00. For 7 = 3 and large enough A'', the crossover 
temperature is found to be Tco ~ ^InA'', with a prefactor A proportional to the mean degree. For 
2 < 7 < 3, we obtain T^o ~ {k)N'^ , with an exponent z that decreases as 7 increases. This exponent 
is found to be lower than predicted by earlier calculations. 



PACS numbers: 64.60.Cn, 05.50.-Fq, 89.75.Hc,84.35.-hi 
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Complex networks describe several kinds of natural 
and artificial systems (social, biological, technological, 
economic), and are currently employed as models to 
study various processes taking place in real-life systems 
[1-3]. In last years, new models of complex networks 
have been introduced, motivated by empirical data in 
different fields. Thus, the so-called small-world [4] and 
scale-free (SF) networks [5] incorporate various aspects 
of real systems. These complex networks provide us with 
the underlying topological structure to study processes 
such as spread of infections [6], signal propagation [1,7], 
and cooperative phenomena [8-11]. 

In a SF network the degree distribution, -P(fc), where k 
is the number of links connected to a node, has a power- 
law decay P{k) ^ k^'' . This kind of networks have been 
found in particular in social systems [12], in protein in- 
teraction networks [13], in the internet [14], and in the 
world-wide web [15]. In both natural and artificial net- 
works, the exponent 7 controlling the degree distribution 
is usually in the range 2 < 7 < 3 [3,16]. 

Cooperative phenomena in complex networks are ex- 
pected to display unusual characteristics, associated to 
the pecuhar topology of these systems. In this context, 
the Ising model on SF networks has been studied with 
several theoretical techniques [9,17,18], and its critical 
behavior was found to be dependent on the exponent 7. 
In particular, when (fc^) is finite, there appears a fer- 
romagnetic (FM) to paramagnetic (PM) transition at a 
finite temperature T^- On the contrary, when (fc^) di- 
verges (as happens for 7 < 3), the system remains in 
its ordered FM phase at any temperature, and no phase 
transition occurs in the thermodynamic limit. 

Here we investigate the FM-PM transition for the Ising 
model in scale-free networks with various values of the ex- 
ponent 7. We employ Monte Carlo (MC) simulations to 
obtain the transition temperature, and compare it with 
that predicted in earlier calculations. Our results con- 



firm those of analytical calculations for 7 > 3, and are 
used to check the precision of those obtained earlier with 
approximate methods for 7 < 3. 

Our networks are defined, apart from 7, by the max- 
imum and minimum degrees, denoted fccut and /cqj re- 
spectively. Thus, the number of sites with degree k 
is given by Nk = (kcut/k)^ for ko < k < fccut, and 
Nk = otherwise. This gives Nk^^^ = 1, and a sys- 
tem size N = X^feo"* ^k, which for large fccut scales as 
N ~ fc^jjj. Once {Nk} is defined [or the corresponding 
probability density P{k) = N^/N], one has a total num- 
ber of ends of links (total degree) K = X^feo"' ^-^fc- Then 
we ascribe a degree to each node according to {Nk}, and 
finally connect at random ends of links (giving a total of 
L = K/2 connections), with the conditions: (i) no two 
nodes can have more than one bond connecting them, 
and (ii) no node can be connected by a link to itself. 
We have checked that networks generated in this way 
are uncorrelated, in the sense that the joint probability 
P(fc, k') for degrees of nearest neighbors fulfills the rela- 
tion P{k,k') ^ kk'P{k)P{k')/{k)'^ [3]. 

On these scale-free networks, we consider the Hamil- 
tonian 



H 



JijSiSj , 

i<j 



(1) 



where Jij = J(> 0) if nodes i and j are connected, and 
Jij = otherwise. S'i = ±1 (j = 1, A^) are Ising spin 
variables. Sampling of the configuration space has been 
carried out by the Metropolis local update algorithm [19]. 
This allows us to study the temperature dependence of 
the magnetization, and in particular the transition from 
a FM to a PM regime as T is increased. Depending on 
the value of the exponent 7, this transition: (i) can oc- 
cur at a well-defined temperature in the thermodynamic 
limit (A'' — > 00), or (ii) may display a FM-PM crossover 
temperature shifting with system size and diverging to 



1 



infinity as ^ oo. For the sake of clarity we will em- 
ploy a different notation for the temperature of phase 
change in both cases. In the first case, we will call it Tc 
as for thermodynamic phase transitions (Tg < oo). In 
the second case, it will be called crossover temperature, 
Tco, to emphasize its size dependence. 
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FIG. 1. Fourth-order Binder's cumulant Un as a function 
of temperature for scale-free networks with 7 = 5 and ko = 3. 
Symbols represent different system sizes: squares, fccut = 15; 
circles, fecut = 17; triangles, fccut = 19, and diamonds, fccut = 
22. These values of fccut correspond to system sizes A'' ranging 
from 4300 to 29300. 

When a FM-PM transition occurs in the thermody- 
namic limit, the transition temperature has been de- 
termined by using Binder's fourth-order cumulant [19] 



Un{T) 



(2) 



where the magnetization M of a given spin configuration 
{Si} is M = Si/N. The average values in Eq. (2) 

are taken over different network realizations and different 
spin configurations for a given network at temperature T. 
In this case, the transition temperature is obtained from 
the unique crossing point for several sizes A''. 

In the second case, the size-dependent crossover tem- 
perature Tco{N) has been determined from the maximum 
of the magnetization fluctuations, (AM)^ = {M'^)n — 
(M)|f, as a function of temperature. We have checked 
that the crossover temperatures obtained by using this 
criterion agree within error bars with those derived from 
the maximum derivative of the heat capacity. In partic- 
ular, both procedures give the same size-dependence of 
Tco in the cases presented below. 

The largest networks considered here included about 
5 X 10^ sites. Such network sizes are required in partic- 
ular to determine the power law characterizing the size 
dependence of the crossover temperature Tco for scale- 
free networks with 2 < 7 < 3. On the contrary, for the 



cases in which a phase transition exists in the thermo- 
dynamic limit (7 > 3), smaller sizes are necessary (see 
below). The results presented below were obtained by 
averaging in each case over 1000 networks, except for the 
largest system sizes, for which 400 network realizations 
were considered. 

Case 7 > 3. — For scale-free networks with 7 > 3, 
the average value (fc^) converges to a finite value as 
fccut — ^ 00. In this case, analytical calculations [9,17] 
predict a well-defined FM-PM transition temperature Tc 
given by 



Tc 2 \{k^)-2{k) 



(3) 



We have calculated Binder's cumulant Un for several val- 
ues of 7 and different network sizes. As an example, in 
Fig. 1 we present Un as a function of temperature for 
7 = 5, fco = 3, and various values of fccut- Tlic transi- 
tion temperature is obtained from the crossing point for 
different system sizes (or cutoffs fccut)- 

The same procedure has been repeated for other val- 
ues of ko. The resulting values of Tc are presented in 
Fig. 2 (open symbols), along with the transition tem- 
perature predicted by Eq. (3) (solid line). In fact, this 
line was obtained by joining values of Tc derived from 
Eq. (3) for integer values of fco in the limit fccut 00. 
The MC results agree; within error bars with the transi- 
tion temperature given by Eq. (3). For comparison, we 
also present the critical temperature obtained in a sim- 
ple mean-field approach [9]: T^-^ = {k'^)/{k), which is 
displayed as a dashed line. 
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FIG. 2. Transition temperature Tc for scale-free networks 
witii 7 = 5 as a function of the minimum degree fco- Sym- 
bols represent results of MC simulations, as obtained from 
Binder's cumulant. Error bars are less than the symbol size. 
The solid line was plotted by connecting points obtained from 
Eq. (3) for integer values of fco and fccut — > 00. The dashed 
line shows the mean- field result for Tc. 
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The critical temperature Tc obtained from Eq. (3) for 
7 = 5 and > 3 can be fitted linearly with good preci- 
sion as Tc = akfi + b, with the parameters a = 1.50 and 
b = —1.72. The value of a can in fact be estimated by 
approximating by integrals the sums giving the average 
values (k) and (fc^) in Eq. (3). This approach gives for 
large fco, Tc « |fco. We note that for ko = 2 and 3, the 
actual values of Tc deviate slightly from such a linear fit. 
For fco = 1, Eq- (3) is not defined, since (fc^)-2(fc) < 0. In 
this case, the simulated networks consist of many differ- 
ent (not connected) components, and Binder's cumulant 
docs not give a crossing for different system sizes. 

Case 7 = 3. — The transition temperature given by 
Eq. (3) increases for increasing 7 and eventually diverges 
for 7 ^ 3, as a consequence of the divergence of (fc^). 
For 7 = 3, analytical calculations [17,18] predict a FM- 
PM crossover at a size-dependent temperature Tco, which 
scales as log N. Such a logarithmic increase of Tco with 
N has been also obtained by Aleksicjuk et al. [20,21] from 
Monte Carlo simulations of the Ising model in Barabasi- 
Albert growing networks. We note that these networks 
have 7 = 3, but display correlations between degrees of 
adjacent nodes [5]. 
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FIG. 3. Crossover temperature Tco for scale-free networks 
with 7 = 3 and fco = 3, as a function of the system size 
N, presented in a logarithmic plot. Symbols indicate results 
derived from MC simulations, with error bars less than the 
symbol size. The dashed line is a least-square fit to the data 
points for N > 500. 

Our results for Tco in the case 7 = 3 and fco = 3 
are shown in Fig. 3 as a function of the system size 
iV in a logarithmic plot. We indeed find a logarith- 
mic dependence of Tco on N, as in earlier works. In 
Fig. 3, symbols represent simulation results and the 
dashed line is a least-square fit to the data points with 
N > 500 (smaller sizes give Tco values that deviate from 
the asymptotic trend). Thus, our results indicate a de- 
pendence: Tc/ J = A\n.N + B, with constants A = 0.83 
and B = —1.28. We have repeated the MC simulations 



for fco = 5 and 9, and obtained the same logarithmic 
dependence as for fco = 3. The prefactor A increases lin- 
early with fco, and in fact we found: A/fco = 0.28 ± 0.01. 
For Barabasi- Albert networks with fco = 5, Aleksiejuk et 
al. [20] found from a fit similar to ours A = 2.6, which 
means A/ko = 0.52. 

For uncorrelated scale-free networks with 7 = 3, Doro- 
govtsev et al. [17] found Tco/. J ~ j{k)\nN. For fco = 3 
and fccut — > 00, one has (fc) = 5.125, and thus their cal- 
culations predict Tco/ J ~ 1.28 IniV, with a prefactor A 
on the order of unity, as that obtained here. Mean-field 
calculations for 7 = 3 give Tc = ^JfcolnTV [9,21], which 
translates into a ratio ^/fco = 0.5, somewhat larger than 
that found from our MC simulations. 
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FIG. 4. Crossover temperature Tco for scale-free networks 
with 2 < 7 < 3, as a function of the system size A'', in a 
log-log plot. Prom top to bottom: 7 = 2, 2.2, 2.4, 2.7, and 
3. Error bars are less than the symbol size. Dashed lines are 
guides to the eye. 

Case 2 < 7 < 3. — For scale- free networks with 7 < 3, 
analytical calculations [17,18] predict a size-dependent 
crossover temperature Tea scaling as ^ J{k)N^, with an 
exponent z dependent of the parameter 7. In Fig. 4 we 
show the temperature Tco as a function of the system 
size N for several values of 7 in a log-log plot, as derived 
from our MC simulations (for fco = 3). The exponent 
7 decreases from top to bottom: 7 = 2, 2.2, 2.4, 2.7, 
and 3. For a given system size, Tco decreases as 7 in- 
creases, as a consequence of the reduction in (fc^). For a 
given 7 < 3 and large-enough networks, log Tco displays 
a linear dependence on log A'', as expected for a crossover 
temperature Tea diverging as a power of the system size 
N. This linear dependence is obtained for system sizes 

^ A'^o, Nq increasing with 7 and eventually diverging 
for 7^3. This means that the present MC procedure 
cannot be applied to obtain accurately the exponent z 
close to 7 = 3, unless one employes much larger system 
sizes. However, for 7 < 2.8, z can be found with enough 
precision for the system sizes considered here. 
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Thus, we have derived the exponent z from our sim- 
ulation results for 7 < 2.8 by obtaining the slope of 
log[Tco/( J(fc))] vs log A'' for large N. Our results are 
shown in Fig. 5 (open symbols) as a function of 7. The 
solid line represents the analytical prediction [17,18]: 



Our results agree with the analytical calculations in that 
the ratio Tco/{J{k)) diverges with TV as a power law, with 
exponent z increasing for decreasing 7. However, our MC 
simulations give values of z lower than Eq. (4) . In par- 
ticular, for 7 — > 2 Eq. (4) gives z = 1, and our numerical 
procedure yields z = 0.43 ± 0.02. 
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FIG. 5. The exponent z giving the power-law dependence 
of the crossover temperature Tco is plotted vs the parameter 
7 for scale-free networks. Open symbols are results derived 

from Monte Carlo simulations. The solid line corresponds to 
the analytical prediction given by Eq. (4) . The dashed line 
represents the dependence z = (3 — 7)77. 

The exponent z in Eq. (4) is related to the divergence of 
(fc^) in networks with ^ P{k)dk ^ l/N, which means 
that fccut ^ A^i/('>'^i) [17,18]. However, for the networks 
considered here, for practical computational reasons we 
have defined a sharp cutoff fccut, which gives a depen- 
dence ^ Ni3-f)h^ or fccut N^^^- Therefore, as- 
suming that Tco diverges with the same exponent as (fc^), 
we have z = (3 — 7)77. This dependence of z on 7 is plot- 
ted in Fig. 5 as a dashed line. Our MC results follow this 
line, but separate from it at 7 = 2. It is not clear the rea- 
son for this discrepancy. One can argue, however, that 
in the limit 7 — > 2 the average value (fc) diverges log- 
arithmically, and thus the convergence of the exponent 
for the ratio (fc^)7(fc) should be very slow, in spite of the 
apparent convergence shown in Fig. 4. This point may 
require further theoretical consideration. 

In summary, we have studied the FM-PM transition 
for the Ising model in uncorrelated scale-free networks, 
by means of Monte Carlo simulations. For 7 > 3 our 



results for the temperature transition fully agree with 
earlier analytical calculations, confirming the appearance 
of a well-defined transition in the thermodynamic limit. 
For 7 < 3 we find a crossover temperature which in- 
creases with system size. In particular, for 7 = 3 such 
an increase is found to be Tco ~ 0.28fcolnA'', whereas 
for 7 < 3 we obtained Tco ~ J{k)N^, with an exponent 
z lower than predicted by earlier analytical calculations. 
We finally note that some care should be taken when 
comparing analytical results for scale-free netwoks with 
those derived from simulations, since the cutoff defini- 
tion for the actually simulated networks may appreciably 
change the results in some cases. 

The author benefited from useful discussions with E. 
Chacon and E. Velasco. This work was supported by 
CICYT (Spain) under Contract No. BFM2000-1318. 
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